3 



Classical nucleation theory from a dynamical approach to nucleation 

James F. Lutsko 1 and Miguel A. Duran-Olivencia 2 

^ Center for Nonlinear Phenomena and Complex Systems, Code Postal 231, 
Universite Libre de Bruxelles, Blvd. du Triomphe, 1050 Brussels, 
Belgium^}. 

2 ) Laboratorio de Estudios Cristalogrdficos. Instituto Andaluz de Ciencias 
de la Tierra. CSIC-UGR, Avda. de las Palmeras, 18100 Granada, 
Spain_l 

(Dated: 17 August 2012) 

»- Eh ' It is shown that diffusion-limited classical nucleation theory (CNT) can be recovered 

as a simple limit of the recently proposed dynamical theory of nucleation based on 
fluctuating hydrodynamics (Lutsko, JCP 136, 034509 (2012)). The same framework 
is also used to construct a more realistic theory in which clusters have finite interfa- 
cial width. When applied to the dilute solution/dense solution transition in globular 
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that have received considerable attention in the recent literature on nucleation: the 
importance sub-dominant corrections to the capillary model for the free energy and 
of the correct choice of the reaction coordinate. 
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I. INTRODUCTION 



The starting point for any discussion of nucleation is the set of ideas collectively known 
as Classical Nucleation Theory (CNT)(see, e.g. Ref. Q). The aim of CNT is to describe the 
growth and dissipation of clusters of a new phase forming within a bath of the mother phase. 
(This paper will primarily address the problem of homogeneous nucleation although the ideas 
discussed here, and indeed those of CNT, are easily extended to the case of heterogeneous 
nucleation.) In CNT, the central quantity is the concentration, cjv, of clusters of the new 
phase of a given size, N. (Note that since we will limit the present discussion to phase 
transitions involving a single species, the number of molecules in a cluster and the mass of 
the cluster are equivalent and the terms will be used interchangeably along with the generic 
term "size".) The basic idea underlying CNT is to write equations for the time evolution 
of the cluster concentrations which are analogous in form to chemical rate equations (the 
Becker-Doring equation^) . The physics of the problem enters via the rate coefficients 
which depend on the properties of the cluster (e.g. its size) as well as on those of the 
mother phase (e.g. its pressure). These equations are then used to extract a nucleation 
rate under the particular assumption of a quasi-stationary distribution in which there is a 
constant concentration of monomers and a steady production of clusters up to the of the 
critical cluster: critical clusters are removed and new monomers added thus allowing for the 
establishment of a steady state. This process is necessary since the system is, by definition, 
not in an equilibrium state and so the concept of an equilibrium distribution is not available. 

One important assumption that is made in CNT, even at this level, is that the concept of 
a "cluster" is well defined. Of course, one can always devise a definition of what constitutes 
a cluster. For example, in their simulation studies of nucleation, ten Wolde and Frenkel 
have defined a liquid cluster forming in a vapor to be the collection of molecules for which 
the number of molecules within a given difference exceeds some threshold - i.e., that the 
local density in the neighborhood of a molecule is above some valued The need for such 
a definition highlights the underlying fact that clusters are not physically distinct objects. 
Indeed, both simulation^ and Density Functional Theory calculations^- show, e.g., that a 
liquid cluster forming in a vapor has a broad and diffuse interface with the surrounding 
vapor. In fact, for small clusters, all molecules in the cluster could be classified as being 
part of the interface with no bulk material present at all. It is only when clusters are large, 
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specifically when the radius of the cluster is large compared to the width of the interface, that 
the distinction between molecules being "within" a cluster and "outside" a cluster makes 
sense. However, this is not a particular problem in the context of CNT as the assumption 
of large clusters is usually made at some point in any case. 

Another subtlety arises in the description of small clusters. Let us consider for a moment 
an under-saturated solution for which there is no difficulty in assuming an equilibrium 
distribution of cluster sizes. In CNT it is usually assumed that the probability density for 
observing a cluster of size N is proportional to the Boltzmann factor P (N) ~ exp (—/3Q (N)) 
where = 1/ksT, ks is Boltzmann's constant, T is the temperature and Q(N) is the 
free energy (e.g. in the capillary approximation) of a cluster of size (mass) N. On the 
other hand, in the capillary approximation, the free energy is often expressed in terms of 
the radius of a cluster, Q(R) = Q(N(R)) and one could equally well suppose that the 
distribution of cluster radii is P (R) ~ exp (— f3Q, (R)\ However, these two expressions 
cannot both be true because, strictly speaking, there is the exact relation P (N) dN = 
P (R) dR from which one would deduce, e.g., P (R) = P (N (R)) j^or, more suggestively, 
P (R) ~ exp (^—f3 (vt(R) + fcgT In dN Jj^ ^- Alternatively, one could take P (R) as being 
more fundamental and deduce P (N) from it with a corresponding shift of the effective 
free energy. Which of these is correct - or are they both wrong? In the CNT limit, these 
distinctions are unimportant since In dN J^ ~ In R and so the correction to the free energy 
is logarithmic and, hence, sub-dominant compared to the volume (R D ) and surface (R D ~ r ) 
contributions that determine Q (R) in D dimensions. However, in the description of small 
clusters, even at low supersaturations, these distinctions become important. One way to 
resolve this ambiguity would be to formulate a Fokker-Planck equation for the probability 
density and then to use this to determine the stationary (e.g. equilibrium) state. This is, in 
essence, the approach followed here. 

Recently, there has been considerable interest in the possibility of extending CNT in two 
different but related directions. The first, discussed in detail by Prestipino et al^ concerns 
the correction of the capillary model for the free energy used in CNT whereas the second, 
discussed e.g. by Lechner et al^, concerns the choice of the reaction coordinate used to 
characterize the nucleation pathway. In both cases, of course, the quantity of most direct 
interest is the nucleation rate. For example, in the former work it was pointed out that 
sub-dominant terms in the expansion of the free energy as a function of cluster size make 
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non-negligible contributions to the nucleation rate for small clusters. The goal of the present 
paper is to demonstrate a similar analysis which has the advantage that the entire nucleation 
processes can be described self-consistently. This is important because it turns out, as 
indicated above, that the choice of reaction coordinate also involves sub-dominant corrections 
to the nucleation rate so that the two avenues of investigation mentioned above are actually 
seen to come together. 

The present analysis is based on a recently developed description of nucleation based on 
fluctuating hydrodynamics^ - — . Its range of validity is restricted to systems governed by a 
diffusive dynamics typical of colloids or of macromolecules in solution. The new formula- 
tion was motivated by difficulties in finding a consistent way to use the tools of DFT (an 
equilibrium theory) to describe the process of nucleation (a nonequilibrium process). In 
this approach, attention is focused on the formation of a single cluster and the fundamen- 
tal quantity is the spatial density distribution that describes the cluster. One feature of 
the theory is that the density distribution can be parametrized in terms of a few physical 
quantities such as some measure of the size of a cluster, the width of the interface, etc. and 
a dynamical description of the evolution of these quantities results. When the cluster is 



i CNT and, indeed, one 
and below. (Subject 
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parametrized by a single quantity, its "size", one makes contact wit 
can reproduce CNT in the appropriate limits as described in Ref. 
to the same approximations such including that clusters are spherical). In this paper, the 
goal is to extend this approach to allow for an extension CNT that includes the finite width 
of the cluster interface in a self-consistent manner. It proves relatively straightforward to 
extract a nucleation rate from the formalism, essentially by following the development from 
CNT, and a comparison between the predictions of the extended theory and of CNT can be 
made. The ambiguity discussed above concerning the equilibrium distribution is resolved 
and it is in fact found that neither the size nor the radius are the most natural reaction 
coordinate. 

In the next Section, the elements of Classical Nucleation Theory are reviewed with a 
particular focus on the assumptions that go into it. The formulation derived from fluctuating 
hydrodynamics is then described. It is noted that the formulation possesses a type of 
covariance with respect to the choice of reaction coordinate so that the ambiguities discussed 
above are resolved. An expression for the nucleation rate is also derived. Section III shows 
how the capillary model for a cluster can be used in conjunction with the dynamical theory 



to reproduce familiar results from CNT such as the expression for the monomer attachment 
rate, the nucleation rate and the rate of growth of super-critical clusters. A model for 
clusters that allows for a finite interfacial width is described in Section IV. Section V presents 
comparisons between the models as well as tests of the assumptions underlying them. Our 
results are summarized in the final Section. 



II. THEORY 

A. Classical Nucleation Theory 

In CNT, one begins by assuming that clusters of a given size can grow by the addition 
of monomers, or can dissipate by spontaneously giving up a monomer to the surrounding 
solution. Processes involving the interaction of clusters which are both larger than monomers 
are ignored. Then, the concentration of clusters of size N, cjv, is determined by a set of rate 
equations of the form 

= Jn-iCn-i — 9nCn + 9n+iCn+i — JnCn (1) 

where /jy is the monomer attachment rate for a cluster of size N and g jy is the monomer 
detachment rate for a cluster of size N. This is known as the Becker-Doring modet^. For 
condensation of a droplet from a diffuse solution, the monomer attachment rate is clearly 
going to be proportional to the rate at which monomers impinge on the cluster and so can 
be written as fx = 1n]n (4vri?^) where Rn is the radius of a cluster of size N, is the 
rate per unit area at which molecules collide with the cluster and 7^ is a phenomenological 
factor describing the probability that a colliding molecule actually sticks to the cluster. The 
rate at which molecules in solution collide with the cluster can be determined by solving the 
diffusion equation with adsorbing boundary conditions with the result that in the long time 
(quasi-static approximation) = Dci/R^ where D is the monomer diffusion constant, 
so that /at = 7tv (^DRn) c\—. The monomer detachment rate is much harder to estimate 
since it depends on the same physical details that are accounted for phenomenologically in 
the sticking constant (i.e. details concerning intermolecular interactions in the condensed 
phase). Kashchiev^ gives a general argument for a relation between the monomer attachment 
frequency and the monomer detachment frequency for an equilibrium (i.e. under-saturated) 
system: namely, that in this case detailed balance demands that J'n-ic^i = 9nc$ so 



that, if we can assume a Boltzmann distribution oc exp (— /3AQ^) then it follows that 
Qn = /at-i exp ((3AQ N — (3AQ N _i). Note that this argument only holds for the under- 
saturated, equilibrium solution and that employing it for the supersaturated, nonequilibrium 
solution is a further assumption. 

In the limit of large ( N >> 1) clusters, Eq.([T]) can be approximated by treating (t) — > 
C (N; t) as a continuous function and expanding to get 

(N, t) _ f c _ 9f { N)C { N;t) + l^f { N)C(N;t) + \ _ c 



dt V dN 2 <9iV 2 

(2) 

where a formal ordering parameter, e, has been introduced: it will, at the end of the calcu- 
lation, be set to one. Simplifying gives the Tunitskii equation^, 

= m ( e {9 (A ° " ; m c {N] t] + V^n (/ (iV) + 9 m c {N] t] ) + (e3) 

(3) 

once we set e = 1 and dropping the third and higher order contributions. Similarly, expand- 
ing the approximation for the detachment rate gives 

g (N) = f (N — 1) exp (/3AQ (N) - f3Att (N - 1)) (4) 

. /W - £ ^ + £/(JV) M) +0(£ *) 

and combing with the Tunitskii equation results in 

iC (n, t) _ , a ff Bf m + f {N) ajA^N)^ c (JV; t) + a ; (JV) c (JV; A + 



dt ' dN VV <9iV J v ' <9iV J K ' ' dN 

-^{fW eJ ^C(N- A+ f ( N)^C(N,t ) ) + 0( < ?) ' (5) 

which is the well-known Zeldovich equation^. 

Although these manipulations are formally correct, there is in fact an additional assump- 
tion being made. Although the derivative ~ can be considered to be small for 
large N, it is not clear that the same is true of the free energy since in fact AQ (iV) — 
/3AQ (N — 1) ~ 0(1). Taking only the first term in the expansion of the exponential in 
Eq.pj is therefore only justified if it can be assumed that 913 <C 1: in other words, in 
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the vicinity of the critical cluster . It is therefore clear that the Zeldovich equation is, strictly 
speaking, only applicable for large N and for TV near the critical size. 

Noting that the probability to observe a cluster of size N, P (N) = ^ n ^n(N') wnere n-(N) 
is the number of clusters of size N, it is easily seen that C (N; t) = C to t (t) P (N; t) where 
Ctot (t) = J2n & (N'i t) is the total number of molecules per unit volume, i.e. the average 
concentration. So, the Zeldovich equation implies a Fokker-Planck equation for P(N;t) 

9 ( SW 9lim W P (N; t) + fm±P (N; t)) - P (N; t) ^ °" "» 



dt <9iV V dN dN ) y ' dt 

(6) 

where the source term vanishes if the total number of molecules is constant. 



B. Generalization based on fluctuating hydrodynamics 

In general, nucleation is a fluctuation- driven phenomena limited by the rate of trans- 
port of mass, heat, etc. within the system. In fluid systems, it is therefore naturally 
described by Landau's fluctuating hydrodynamics which includes both transport and ther- 



mal fluctuations^ 1 ^, 
described in Refs. 



10 



he description of nucleation via fluctuating hydrodynamics has been 
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A particularly simple description is possible when the system is 
diffusion limited and over-damped, as is an appropriate description for colloids and macro- 
molecules in solution. The dynamics can be reduced to an equation for the time-dependent 
spatial density and a reduced description is possible when the density field is parametrized 
in the form 

p{r;t)^p{r;X x {t),X 2 {t)...X N {t)) (7) 

so that the time dependence occurs via a set of N parameters which could include, e.g., 
the size of the cluster, the width of the interface, the density at the center of the cluster, 
etc. Here, we will be interested in the case that there is a single parameter. Then, it was 
shown^— that the evolution of the order parameter may, in the case of an infinite system, 
be approximated by a stochastic differential equation of the form 

f^-vwf-^w^f+vwwfw (8) 

where D is the diffusion constant, Q(X) = F (x) — pN (X), is the grand potential for 
chemical potential p and where N (X) is the total number of molecules. The last term 



on the right is proportional to £ (t) which is a delta-function correlated fluctuating force 
with mean zero and variance equal to one. The quantity g (X) that determines the kinetic 
coefficient of the SDE as well as the amplitude of the noise is given by 

r 1 fdm(r;X)\ 2 , 



where the cumulative mass is 



pr 

m(r;X) = 4ir p (r; X) r' 2 dr' . (10) 
Jo 



'0 

Since the amplitude of the noise in the SDE is state-dependent, it is important to note 
that the equation is given according to the Ito interpretation. The function g(X) will be 
called "the metric" as it can be used to define a meaningful length in the general theory^ , 
although it will not be used for that purpose here. 

1. Fokker-Planck equation 

It is useful to also consider the equivalent Fokker-Planck equation which is^ 

Note that this can be written as 

= D JL U (x) s {B W - '"^ W) + m 4.) p { x, t) (i2) 

dt dX y y K J OX y K 1 dX J K ' J K J 

so that comparison with Eq.Q shows that these are formally the same with the monomer 
attachment frequency replaced by g^ 1 (X) and with the free energy shifted by a term log- 
arithmic in g{X). From the expression for g(X), Eq.flS]), it is evident that this function 
varies as the volume of the cluster so that the shift to the free energy will go like the log of 
the radius of the cluster. 

Let us look for a stationary solution determined, for some constant J s , from 

- D (g-' (X) 9 -^l + g-W (X) ±g~W (X)) P s (X) = J s , (13) 

giving 

P s (X) = Ag 1 ' 2 (X) exp (X))-^- 1 J s g 1 ' 2 (X) exp (-f3Q (X)) / g 1 ' 2 (X') exp (f3Q (X')) dX' 

(14) 



for some constant A. Now, if we note that for some change of variables X — > Y (X) we have 
that 

dX x ' 



9(Y)=g(X(Y))( — ) (15) 



and so 

»y 



P s (Y) = Ag 1 ' 2 (Y) exp {-/3Q (Y)) - D-\j s g l l 2 (Y) exp {-(3Q (Y)) J g 1/2 (Y') exp {(3Q (Y')) dY' 

(16) 



dX 

= AgV 2 (X)—exp(-pn(X)) 

dX f x dX' 

- D~ 1 J s g 1 / 2 (X) — exp (-fiCl (X)) J g 1 / 2 (X') — exp (X')) dY' 

Hence, this solution is completely covariant - there is not ambiguity as to which variable 
is used. If the system is stable (i.e. under-saturated), then it makes sense to seek an 
equilibrium ( J = 0) solution which is 

P eq (X) = Ag 1 / 2 (X) exp (-/3fi (X)) = A exp (-0 (n (X) - h B T\ng (X)\ \ (17) 

where the constant A is fixed by normalization. Note that this has the form of a canonical 
distribution in terms of the shifted free energy. 

2. A Canonical Variable 

In the present, single-variable, case there is always a special variable for which these 
expressions simplify. Given any variable, X, it is defined via 

dY = y/g (X)dX (18) 

and an arbitrary boundary condition that will be taken to be Y (0) = 0. In terms of this 
variable, the Langevin equation becomes 

where Q (Y) = Q (X (Y)). The corresponding Fokker-Planck equation is 

BP(Y,t) 3 f e0a(Y) B\ ? 

a ~ dY\ dY + dY) { '>- { ' 



These are equations that might have been written down on phenomenological grounds but 
they would have been ambiguous since there is no a priori reason to use the mass as the 
independent variable rather than, say, the equimolar radius. This illustrates the way in which 
grounding the theory on a more fundamental description serves to remove such ambiguities. 



3. Nucleation rate 



The standard argument to obtain the nucleation rate has a long history, as discussed by 
Hanggi et aU^, and is based on boundary conditions according to which (a) the distribution 
is stationary and (b) the distribution goes to zero at some point beyond the critical cluster, 
say at X = X + (see also Refs. (l| and 17)). The latter condition represents the physical 
fact that once a cluster is slightly larger than the critical cluster, it will almost certainly 
grow forever and so plays no role in the stochastic part of the process. One therefore 
imagines establishing a steady state by removing clusters once they reach size X = X + and 
simultaneously re-injecting the removed material in the form of monomers. Here, since we 
are only interested in the case of one-dimensional barrier-crossing, this rate is easily evaluated 



from the exact solution for the steady-state distribution 
approximations are necessary as described in Refs. Il6|, 



or multivariate problems, further 



18, and 



19|) 



First, the nucleation rate, J, is the rate of production of super-critical nuclei per unit 
volume and is therefore 



J 



d 



-n (N; t) dN 



d 



N c (t) 
V 



P (N; t) dN 



(21) 



dt J n* V dt j n 

where N c (t) = J °° n (N; t) dN is the total number of clusters. Making use of the Fokker 
Planck equation, this becomes 

1 dN c (t) 



J 



V dt 

N c (t) 



P (N, t) dN 



(22) 



D 



V 



N*+l 



g- 1 (AO 



df3Q (N) 
dN 



P (N, t) + g-W (N) (AT) P (N, t)j 



N=N* 



In the artificially imposed stationary state, the first term on the right does not contribute 
leaving 

df3fl (iV) 



J = -Df (V 1 (AO ^§^Ps (N) + g-W (N) ±g~W (N) P s (JV)) 



(23) 



N=N* 
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Since we remove clusters of size X + as they form, we need a steady-state distribution that 
satisfies P (X + ) = 0. When this is imposed, the general steady-state distribution, Eq.f JT^|) . 
becomes 

P s (X) = D^Jsg 1 ' 2 (X) exp (-00 (X)) C + g 1 ' 2 (X') exp (X')) dX' (24) 

Jx 

Substitution into the expression for the nucleation rate gives 

N 

J = yJs (25) 

The remainder of the development concerns the relation between the imposed stationary 
flux, J s , and real properties of the system. Zeldovich, as described by Kashchievi requires 
that the concentration of monomers be fixed at the "equilibrium" value 

Poo = yPs (Xi) (26) 

where X\ is the value of the parameter X that fixes the number of molecules in the cluster 
to be 1. (The notation chosen here is motivated by the fact that one expects that Pmfty is 
the density far from a cluster.) This then gives 

Jcnt = v~ (27) 

gW (X,) exp (-pn (X 1 )) f£ gW (X') exp (/3Q (X')) dX> 

Dpoo 



/*+ g (X') exp (A/K) (X') - \ In jjfj) dX> 

The second equality gives the result written in terms of the adjusted free energy. A saddle- 
point evaluation for the choice of number of molecules as the variable gives 



J CNT ~ D Poo g~ l (AQ J^-y/\A/3n" (AUIexp (-A/3fi„) 

V Z7T 



(28) 



where the adjusted critical size is determined from 



and where 



A(3Q^ = A/3Q (N„) - ~ In ^4 ( 30 ) 

1 9{ A i) 

The problem with this expression is that it is not covariant and the result will depend on 
the choice of X. This is due to the fact that the condition, fixing the fraction of monomers, 
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is itself not covariant since only the combination P s (X) dX is invariant under a change 
of variables. We stress that within the context of CNT this fact is irrelevant since the 
assumption of large N means that the lack of covariance is due to sub-dominant terms in 
the distribution. However, within the context of the general theory, this lack of covariance 
indicates unnecessary ambiguity since different results will be obtained for different choices 
of the variable X appearing in the boundary condition, Eq.( 126]) . 

The first step to solving this problem is to note that, as defined here, P s (X) is normalized 
so that we can immediately fix the unknown coefficient giving 

p , x) = 9 1/2 (X) exp (-gft (X)) {X') exp (gft (X')) dX> 

f X+ dXgVHX)e Xp (-{3n(X))fX + g^(X>)e W m(X'))dX' 

and 

j — t (32) 

V /*+ dX g 1 / 2 (X) exp (-/3ft (X)) /*+ g 1 / 2 (X') exp (/3ft (X')) dX' 

This still leaves the question of determining the overall concentration of clusters. We do this 
following the general idea used in CNT but being careful to preserve covariance by imposing 
a condition on the total number of molecules, N" (Nq) in clusters up to size No, 

/-No rNo 

M (N ) = Nn (X) dN — N c NP S (X) dN (33) 
Jo Jo 

so that 

N _ y , N ) Io X+ dX 9 1/2 (X) exp (-/3ft (X)) + gV* (X') exp (/3ft (XQ) dX' 

J X{No) dXN(X)g 1 / 2 (X)exp(-^n(X))J^ + g 1 / 2 (X')exp(pn(X'))dX' 

(34) 

giving the nucleation rate as 

J(No) = 

j* {No) dX N (X) g 1 / 2 (X) exp (-/3ft (X)) j* {N+) g 1 / 2 (X') exp (/3ft (X')) dX' 

(35) 

where the average density is p (X ) = A/ "^ v °' ) and where we have indicated that the boundary 
condition makes the nucleation rate a function of Nq. If N is chosen to be large, say the size 
of the critical cluster X*, then Af (N*) is the total population in the initial (i.e. precritical) 
state and this is essentially the standard expression for the reaction rate derived from one- 
dimensional barrier crossing (the "flux over population" expression, see, e.g., Hanggi et at^) 
except that an extra factor of X occurs in the integral in the denominator, which is meant to 
count the population in the initial state, because each cluster contains many molecules. If X 
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is chosen to be on the order of one, then we are essentially counting the number of monomers. 
Importantly, if most of the material exists in the form of small clusters (monomers, dimers, 
etc.) then these estimates will be the same. If they are not the same, then a significant 
amount of mass is present in the form of clusters and the implicit assumption (made here 
as well as in CNT) that clusters do not interact is probably invalid. One of the goals below 
will therefore be to monitor the validity of this assumption. 

Assuming that the bulk of the material is in the form of small clusters and that Nq is 
chosen sufficiently large , an approximation to the exact expression, Eq. (1351) . can be devel- 
oped as described in Appendix |A] Assuming that the free energy and number of molecules 
as functions of the canonical variable Y have the expansions Af3Q(Y) = QoY a + ... and 
N (Y) = uqY 13 + ... for small Y, it turns out that the approximation implies that, for small 
X, the stationary distribution can be approximated by 

P s (X) ~ f^g l/2 (X) exp (-A/3Q (X)) (36) 

and in which case the nucleation rate is approximated by 

1 „ afi ° 



PW (X,) g- 1 (X,) exp(-A/3f2.) (37) 



III. CAPILLARY MODEL: CLASSICAL NUCLEATION THEORY 

The process we will describe is the nucleation of a liquid, with bulk density pi, from a 
vapor with bulk density p v at some temperature T and chemical potential p,. Let u (p) be 
the free energy (grand potential) per unit volume, so that u (p) = / (p) — pp where / (p) 
is the Helmholtz free energy per unit volume. Then, the liquid and vapor densities are 
determined by the imposed chemical potential via u' (pi) = u' (p v ) = and since we choose 
thermodynamic conditions such that the vapor is metastable, u (pi) < u(p v ). To recover 
CNT it is only necessary to use the capillary model for the density 

{Po, r < R 
(38) 
Poo, R<r 

where we take p = pi for the density inside the cluster and p^ = p v for the density outside 
the cluster. The capillary-theory expression for the free energy of the cluster is 

A-jr 

A/3Q (R) = —R 3 Af3u (p ) + 4?ri? 2 7 (39) 
13 



where the second term represents the effect of surface tension. Note that the only parameter 
that is allowed to vary is the radius of the cluster. 

The model for the density gives the cumulative mass density 

t ^ I T^o, r < R 

m(r) = < 6 (40) 

1 f#Vo + f (r 3 -i? 3 ) Poo, R<r 

and the metric 

9 (R) = r -j-4— (^R 2 (Po - Poo)) 2 dr = (P °- P °° )2 4rf (41) 
J R 4vrr^poo Poo 

The excess number of molecules in the cluster is 

4:71 

AN = —R 3 (p - Poo ) (42) 
and the metric in terms of the number of molecules is 

9 (AJV>= 9 (*)(^L)' = — -L^. (43) 
The canonical variable is 



Y = I J (P ° ~ p ^ AnR'W = -W47r (Po ~ Poo) V 2 . (44) 



In comparing the Fokker-Planck equation for the generalized theory, Eq. (ll2[) . to that of 
Zubarev, Eq(jS]), one finds that agreement provided that the logarithmic corrections to the 
free energy are neglected and the effective monomer attachment frequency is identified as 

/ (iV) = g" 1 (AN) = 4ir Poo R (AN) (45) 

which is the usual result for diffusion-limited nucleation in the case that the phenomenolog- 
ical "sticking constant" is equal to one^. The stochastic differential equation for the radius 
is 

dR Poo d(pn-\ng 1 l 2 (R)) , — - fS 

When the cluster is large (i.e. when it is super-critical) the noise becomes unimportant and 
the radius grows as 

dR = D Poo \A^(p )\ R ^ + Q (47) 

dt (p - Poo) 
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which gives the classical result R ~ t 1//2 when the higher order terms are neglected^. In the 
weak liquid limit, Poo <C po, the coefficient of agrees with that given by Lifshitz et at^. 
The (non-covariant) CNT-like nucleation rate, from Eq. (l28|) . is 



Jcnt ~ 2V2^D Poo R (AN,,) v/|A/3ft" (AAQ|exp (-A/30**) (48) 

and in the limit that the logarithmic corrections to the free energy are negligible, this agrees 
with the result from CNT. In the following, we will use as a reference the usual CNT result 
that is obtained by ignoring the logarithmic shift in the free energy, 

Jcnt ~ 2V^D Poo R (AN,) ^/\Af3Q" (AN*)\ exp (- A/3fi™ T ) (49) 

Po y/l 3 (A(3u) 

The approximate nucleation rate based on a condition of fixed mass, Eq (l37[) . is 



J 3p v / 2(4 7 r)^ 13 I n , , i , \ 1 1 / 2 

-j rfu , - — 77- R 7 4 \pn (R*)g (R*)\ (50) 

Jcnt T (—) ^fp^ \Apu\ (p - p^) 

3V2(47r)t p o 1/2 

1 (-J (PO - Poo) 

Note that in all cases, we are assuming that the average concentration of material in the 
(hypothetical) steady state, p, is the same as the background density, Poo . 

IV. EXTENDED MODEL: FINITE CLUSTER WIDTH 

A. The cluster structure and the metric 

A significant short coming of the capillary model is that the width of the cluster's interface 
is zero. A more realistic model will have a finite width which, from various simulations and 
DFT calculations, might be expected to be two or three molecular diameters in width. A 
simple extension of the capillary model to take account of a finite width is the piecewise- 
linear profile, 

Po, r < R — w 

P(r) = I Po - (Po - Poo) R-w<r<R (51) 

Poo, R<r 
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The corresponding cumulative mass distribution for < R < w is 

m (r) = e{R-r)^- (p - Poo ) r 3 (AR - 3r) + 6 (r - R) (p - Poo) R 4 + V (r) Poo (52) 
while for R > w it becomes 



m (r) = (p - Poo) V (r)Q(R-w -r) 



(53) 



+ 6 (r - (R - w)) 9 (R - r) (p - Poo) (r 3 (4fl - 3r) - (fl - w) 4 ) 

3w 



+ (r _ ^ Z. (po _ Poo) (ij* _ (ij _ 

+ V(r) Poo . 



Calculation of the metric is then straightforward with the result that for < R < w 



g{R)=^ 



with 



PO - Poo 

3w 



wR 4 



.(Po ~ Pc 



a 4 In 



a = 1 + 



a — 1 



u>p c 



a 3 H — a 2 H — a H — 



1 1 
— a H — 
3 4 



Poo 



R (PO - Pc 

Note that for small clusters, R <C w, this gives 



9(R) 



Itt ( po - poo V i? 5 



15 \ w / Poo 
For larger radii, R > w, the result is 

■ 

0(A) 



R 



1 _ (PO - Poo) ~ Q 

36poo to 



(54) 
(55) 

(56) 



4vr (po-Poo) 2 {R 3 -(R-w) 3 Y 



9 p c 



2 i? 



(57) 



4tt / p -poo \ / (g - w) b R 
9 V I \ R 2 a 2 U R-w 



((R 



w 



R 3 a 3 )' 



R 2 a? 



In I ^ 



^_ ( P0~ Pc 

9 



-(2 (*-»)» 



Po 
i? 3 a 3 ) iw 



-±i? 2 a 2 (# 2 - (i? - u>) 2 ) + |i?a (i? 3 - (R - wf) + 1 -(R i -(R- w) 4 ) 

One can again define a canonical variable using 

dY 



dR 



and for small clusters one finds 



= 2 / Po -p c 
" 7 1 w 



15p 



i^(l + 0(i?)). 



(58) 



(59) 
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B. Free Energy model 



It would be somewhat inconsistent to use the capillary approximation for the free energy 
given that the assumed profile now has finite width. We therefore consider a simple, but 
more fundamental free energy model based on the squared-gradient approximation, 

il\p) = J L (p (r)) + l -K (Vp (r)) 2 ) dr (60) 

where u (p) = f (p) — pp is the grand potential per unit volume, / (p) is the Helmholtz free 
energy per unit volume for a bulk system with uniform density p which can be determined 
based on a given pair potential using thermodynamic perturbation theory or liquid state 
integral equation methods. The squared-gradient coefficient, K, can be estimated from a 
model interaction potential using the results of Ref. 21] . For the assumed density profile, 
this becomes 

Air q 
Q (R; w)-tt 00 = — (max (R - w, 0)y Au (p ) (61) 
3 

+ 4tt / Aw U - (p - Poo) L r 2 dr 



max(0,i?— w) 



+ \k^- (R 3 - max (R - w, 0) 3 ) ( 



PO - Poo 
W 



W 

2 



The result for R> w, can also be written as 



Q {R; w) - ^ = — (R - wf Aw (p ) (62) 
+ 4n \ 2u w + K ip0 ~P°° )2 )R* 

4tt I 2ZJ lW + K ^ '^ 2 | Rw 



\2 



T ^{P0 — Poo) 

+ 4:71 | u 2 w + K I w 

6w 



2 



where the density moments of the excess free energy per unit volume are 

1 



U n = ~ / (W (x) - U (Poo)) (X ~ Poo)" dx, (63) 



(PO — Poo)" +1 J poo 

thus showing how, for large clusters (in the ^ — > limit), one recovers something like the 
capillary approximation but with a variable width. Minimizing with respect to the width 
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at constant radius and solving as an expansion in the radius (i.e. assuming R 3> w) gives 
(see Appendix [B]) 

A , A ^(Po) -2ZJiw;o \ f 

w - = Wo { 1+ A^^R + --j (64) 

with 



^(Po - Poo) 2 

""^(Bb-^fo)) (65) 
and the free energy becomes 

fi(i2;w)-n 00 = (^ 3 ) Au(p )+{4nR 2 ) (±Au (p ) + - (P ° ~ q Poo) j +0 {R 1 ) (66) 

The higher order terms are a simple illustration of the post-CNT corrections to the free 
energy barrier recently discussed by Prestipino, et al&. At lowest order, the implied capillary- 
like model is 

n (R; w)-n oc = (y r3 ) Alo 0°°) + ( 47ri?2 ) 7 ( 6? ) 



with 



1 = 1 -Au J ( Po )w 0+ Kip °- p - )2 (68) 



Auj (p ) 



2if (p - poo) 2 (u - Au (p )). 



4 (wo - Aw (po)). 

Note that the coefficient 7 depends on the supersaturation implicitly via its dependence on 
the densities. In CNT, it is more common to ignore the state dependence of this coefficient 
and to fix its value to that which gives the correct value of the planar surface tension at 
coexistence. Here, "correct" will be taken to mean that it gives the same value as the full 
model with finite width. Since at coexistence the free energies of the two phases are equal, 
this is 



Icnt = (P c ° ex ~ PZ ex ) V^mr. (69) 

where the superscripts indicate that the densities are those at coexistence. 

One peculiarity of this model is that the value of the width that minimizes the free energy 
undergoes a bifurcation as the radius increases. This is due to the fact that for w < R the 
free energy is 

n (R; w)-n oD = 4vr J* Aco (p - (p - Poo ) r ~ ( ^~ W) ) r 2 dr + ^-KR 3 ( Pi ^) ' 

(70) 
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and it is clear the free energy difference decreases monotonically to zero as the width in- 
creases. 

Expanding in R one has that 

(R; w)-Q eo = %K (^LZPsA 2 i?3 + ^ Aajl> ( } fPo^Poo \ 2 R5 + (71) 
3 \ w J 5! \ w J v/ 

and in terms of the canonical variable 

n { Y-, w )-^ = ^K(e^Y (U^L.),/^] Yi+OIYT) (72) 



3 \ w J \ 2 VPo-Poo/ V 87T 
which can be used to evaluate the approximate nucleation rate 

J ext ~ Z ext \pn" (R*)\ 1/2 exp (-00.) (73) 

with the prefactor 

Z ext = 6 ^K^(cJ - Auj (po)) 1 g~ 1 ' 2 (R*) (74) 

V. RESULTS AND COMPARISONS 

In order to illustrate the theory developed above, we have performed detailed calculations 
for a model globular protein. We assume that the solvent can be approximated, crudely, 
by assuming Brownian dynamics of the (large) solute molecules which also experience an 
effective pair interaction for which we use the ten Wolde-Frenkel interaction potential 

foo, r < a 
u ( ( i \ 6 ( i \*\ ^ ^ 



2 -i V(^) 2 - 1 



- « Tr^T ) , r>a 



with a = 50 which is then cutoff at r c = 2.5a and shifted so that v (r c ) = 0. The temper- 
ature is fixed at /c#T = 0.375e and the equation of state computed using thermodynamic 
perturbation theory. The transition we study is that between the dilute phase and the dense 
protein phase which, in the present simplified picture, is completely analogous to the vapor- 
liquid transition for particles interacting under the given pair potential. Throughout this 
Section, the supersaturation will be defined as the ratio of the density of the vapor phase to 
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FIG. 1. The phase diagram for the model globular protein. The full lines show the dilute/dense 
solution coexistence, the broken lines are the spinodal and the horizontal line shows the temperature 
used in the calculations. 

that of the vapor at coexistence, S = p v /p vc - The gradient coefficient, K, is calculated from 
the pair potential using the approximation given in Ref.— 



n rj roc 

(3K ~ ~^d 5 (3v (d) + — J {2d 2 - 5r 2 ) (3v (r) r 2 dr 



(76) 



where d is the effective hard-sphere diameter for which we use the Barker-Henderson ap- 
proximation. For the temperature used here we get j3K = 1.80322a 5 . 

In the following, we will characterize the size of a cluster by either the total excess number 
of molecules in the cluster, 

AN = J (p(r) - Poo ) dr (77) 

or by its equimolar radius, Re, which is related to AN by 

Air 



AN = ^-R 3 E (p -p c 



(78) 



For the simple capillary model one has that AiV = ^i? 3 (po — Poo) and Re = R- For the 
extended model, a simple calculation gives 



7T 



AN = -— (R A - (max(i? - w, 0)) 4 ) (p - Poo ) 

6W 



(79) 



A. Energy of formation of a cluster 



The dilute-solution/dense-solution phase diagram is shown in Fig. [TJ By definition, 
the coexistence concentrations correspond to saturation S = 1 and at the spinodal the 
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FIG. 2. The free energy as a function of cluster size (excess number of molecules in cluster) at 
5 = 1.025, 1.175 and 2.5 as calculated using the capillary model with fixed 7, Eq. (139p and Eq. (|69p . 
the capillary model with variable 7 (Eq. (|68p ). the extended model with fixed width, Eq. (|6ip and 
Eq. (]65p . and the extended model minimized with respect to the width. 



supersaturation is found to be S = 2.99. We will therefore illustrate the results of the 
models for supersaturations from S = 1.025 to S = 2.5 corresponding, as will be seen, to 
quite large and very small critical clusters, respectively. 

Figure [2] shows the energy of formation of a cluster as a function of its size for different 
values of supersaturation. The figures show the energy as determined using CNT (capillary 
model with the value of 7 = ■Jcnt fixed to give the correct surface tension at coexistence (S = 
1), Eq. (j69~j) ). the same model but with a supersaturation-dependent value of 7 (Eq. (l68l) ). 
the extended model with fixed width and the extended model minimized with respect to the 
width. It is clear in all cases that the capillary model with fixed 7 gives lower free energies 
than when 7 is allowed to vary and that it is also in closer agreement with the finite-width 
model. This is somewhat counter intuitive. On the other hand, the extended model with 
fixed width gives virtually the same results as when the energy is minimized with respect to 
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the width, showing that the simple fixed-width model is adequate. For these reasons, only 
the capillary model with fixed 7 = jcnt and the extended model with fixed width will be 
used below. 



B. The stationary distribution 

The exact stationary distribution is given in Eq. (l24p . Although both the capillary model 
and the extended model depend on a single variable, a "radius" denoted R in both cases, 
the meaning of the parameter is not the same. In order to make a meaningful comparison, 
we therefore give the stationary distribution using the equimolar radius as the independent 
variable and noting that 

P(R E ) = P{R)^r- (80) 

For the capillary model R = Re so that P(Re) = P{R)- For the extended model, we find 
that 

dR E R 3 - (max(R - w, 0)) 3 



dR 3wR% 



(81) 



Figure [3] shows the stationary distributions as calculated in CNT and using the extended 
model for different values of the supersaturation. Also shown are the approximate distribu- 
tion used to evaluate the nucleation rate. It is apparent that the approximate distribution 
has the right shape so that it is indeed the case that the stationary distribution is well 
approximated by the "equilibrium" distribution, g(X)e~ l3An . However, in the case of the 
extended model, the approximate evaluation of the normalization of the distribution is poor 
due to the rapidly changing analytic structure of the free energy as a function of cluster 
radius for small clusters. A surprising result is that the distributions for the capillary model 
and for the extended model differ even for small supersaturation. This is simply a reflection 
of the fact that the differences induced by the two models are most pronounced for small 
clusters, regardless of the supersaturation. 

The approximate forms for the nucleation rates given in Eqs. (1501) and (1731) are not cal- 
culated from the exact stationary distribution but, rather, from the approximation given 
in Eq. flA6p . The figures show that for CNT this approximation is quite good at low su- 
persaturations but is in considerable error for small clusters. The reason for the increasing 
error is illustrated in FigJH which shows the convergence of the "exact" expression, Eq. (l35l) . 
as a function of the domain of integration. For large supersaturations, the convergence is 
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FIG. 3. The stationary distribution for S = 1.025, left panel, and S = 2.5, right panel. Exact 
results, i.e. evaluations of Eq. (|24p are shown for the capillary and extended models in units such 
that <7 = 1. Also shown is the approximation, Eq. (|36j) . and the approximate distribution scaled so 
as to have the same maximum as the exact result. 



slow indicating that larger clusters are contributing significantly to the evaluation of the 
nucleation rate. For large supersaturations, this calls into question both the monomer at- 
tachment picture that underlies CNT and the assumption of non-interacting clusters that is 
tacitly used in both CNT and the dynamical theory. The main conclusion to be drawn from 
this result is that even if one could accurately calculate the free energy and distribution of 
clusters at large supersaturations, it is most likely the case that cluster-cluster interactions 
would invalidate the assumptions necessary to calculate the nucleation rate, both in CNT 
and in the general framework described here. 



C. The nucleation rate 

Table [I] gives the nucleation rate as calculated for the two models. For the lowest super- 
saturation, the extended model gives a much higher nucleation rate. This is almost entirely 
attributable to the fact that the free energy of the critical cluster is about 10k bT lower in 
the extended model than in the capillary model (and, indeed, the log of the ratio of the 
nucleation rates is about 9.7). In all other cases, the nucleation rates of the two models are 
comparable and similar to that given by the classical CNT formula. 
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10 20 30 10 20 30 

N N 

FIG. 4. The convergence of J(N), Eq. (|35j) . as a function of A for the capillary model (left panel) 
and the extended model (right panel) for different values of the supersaturation. Note that in the 
case of the highest supersaturations, the nucleation rate does not reach a stable value before the 
integral extends up to the size of the critical cluster, A*. 

TABLE I. Properties of the capillary and extended cluster models as functions of the Supersatura- 
tion, S. The classical nucleation rate, Jcnt is given followed by the critical equimolar radius, R*e, 
the excess mass, A A*, and excess free energy, A$7*, of the critical cluster. Also given are the nu- 
cleation rates as calculated from the "exact" expression, Eq. (|35p . and the approximate expression, 
Eq. (|37p . J and J app , respectively. Lengths are in units of a and times are in units of a 2 /D. 
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VI. CONCLUSIONS 

In this paper an alternative approach to the semi-phenomenological description of nu- 
cleation has been presented. The theory is represents a simplification of a description of 
nucleation based on fluctuating hydrodynamics. One of the advantages of this approach is 
that there is no ambiguity regarding the choice of reaction coordinate. It also allows for the 
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incorporation of more complex cluster models. When evaluated using the capillary model for 
a cluster, the results of Classical Nucleation Theory were recovered. A more realistic model 
that includes a finite interfacial width was also developed. This gives rise to sub-dominant 
corrections to the free energy of the type described by Prestipino et al^. It was also shown 
that the different choices of reaction coordinate give rise to corrections to the free energy 
that scale logarithmically with the cluster size. When the different models were used to 
evaluate the nucleation rate, it was found that the dominant effect was the difference in the 
free energy of the critical cluster in the two models. Thus, despite the fact that the model 
for the structure of a cluster affects all aspects of the dynamics in our general approach, it 
is nevertheless the commonly assumed, that only the correction to the free en- 

ergy barrier is really important in determining the nucleation rate. This provides post-hoc 
justification for the combination of the framework of CNT with the use, e.g., of DFT to get 
better estimates for the free energy barrier. 

We also were able to explicitly calculate the stationary distribution. It was shown that 
the stationary distribution is well approximated, for small clusters, by the "equilibrium" 
distribution. However, at high supersaturations, it was found that the stationary distribution 
becomes quite broad and that the nucleation rate becomes dependent on the domain of 
integration over the cluster sizes. This indicates that the assumption of cluster growth 
by monomer attachment and detachment is probably invalid as is the assumption of non- 
interacting clusters. It is therefore the case that quantitative prediction of nucleation rates 
is not really possible with the theory discussed here. For the conditions considered here, 
the breakdown occurs when the critical cluster contains between 20 and 100 molecules. A 
cluster of 100 molecules is still very small and our results illustrate, in some detail, the fact 
that Classical Nucleation Theory is remarkably robust given the crude assumptions that go 
into it. 
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Appendix A: Approximate Nucleation Rate 



The nucleation rate is 



J- 1 = p^D- 1 [*' N (X) g 1 / 2 (X) exp (-/3Q (X)) (J* + g 1 ' 2 {X') exp (/3Q (X')) dX'^j dX 

(Al) 



and is the exact expression for the nucleation rate given the assumptions of the theory. In 
order to develop an approximation it is convenient to switch to the canonical variable so 
that the nucleation rate becomes 



J 1 = 



* N (Y) (y)) (J + exp (pn (y')) dY') dY (A2) 



The first step is to assume that the free energy has a maximum at y which is defined by 



dn (y) 



dY 



(A3) 



so that we can evaluate the inner integral by expanding 



/3fi(y) + i/3fi" (y)(y'-y) 2 



(A4) 



giving 



J^-p-^exp /3fi(y 



\ 



2n 



f3Q" (y; 



J N(Y)exp(-[3tt(Y)^dY (A5) 



Before proceeding, note that this is really an approximation to the stationary distribution 
of 



p (y ) ~ jzr 1 exp ( pn (y 



2tt 



exp (-(3Q (y) 



(A6) 



or, in general, 



P (X) ~ JZT 1 exp (/HI (X,)) J 2 ^!^ exp ( -/HI (X) 



(A7) 
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To further simplify, we assume that the free energy has a minimum at Y = and that 
f3Q (Y) = aY a + ... and N (Y) = N Y^ + ... for some values of a, (3 > giving 



or 



j- 1 ~ p^D^exp^n (y 

= p^D~ 1 exp (fit (y 



2tt 



f NqY 13 exp (-aY a ) dY 
Jo 



(At 



\ W'(Y 



27T AT 1 - m f°° y & + 1 - a I V\ A7 

JM —a <* I Z <* exp{—Z)aZ 
« Jo 



\ PW{Y, 



N Q —a « 1 

a \ a 



J ~ Dp 




1 i+£ a 

1 2^ a Q Fp±^y 



1/2 



/3ft" (y) exp -/3ft (y 



(A9) 



It is easy to translate into the original variables to find 



J ^ Dp 




1 i+£ a 
' 2^ a " f^HJ 



/3ft" (R*) g~ l (RJ 1 1/2 exp (-/3ft (i^)) . (A10) 



Appendix B: Expansion of the free energy 



Beginning with the model expression for the free energy, 



4vr , , 3 



n(R;w) -fioo = — (R-wYAuj(po) (Bl) 

w 

4zj iW + A- (Po ~ Poo)2 ) 
to 

o— i 1 r^(P0 — Poo) 2 i 

3 w 



2 



we minimize by setting the derivative with respect to the width, w, equal to zero giving 

- wf Aoo (po) + ( 2uJ - K (P°~P°°) | b 2 - 8uJ lW R + 6uJ 2 w 2 + \k (p - Poo ) 2 

(B2) 
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or, after dividing by R 3 and rearranging, 



= ((-2Aw (p ) + 2w ) w 2 - K (p - Poof) 



(B3) 



w 



+ (4Aw (po)-8uJ 1 )w 2 - 

JrC 



+ (^-2Aw (p ) w 2 + 6lo 2 w 2 + (p - Poo ) 2 ^j ( 



w\ 2 



We now expand in inverse powers of the radius, 



w 



w = w | 1 + ax— + a 2 



Wo\< 

RJ 



(B4) 



and solve order by order in ^ giving 



K (p - p c 



a 2 



and 



" 2 (-Aw (p ) + w ) 
Aw (p ) - 2cJi 



Aw (p ) - w 
1 



6 (Aw (p ) - w 



—2 (11 (Aw (p )) + Aw (p ) (5w - 60wi + 9w 2 ) - 9u m ~ 



Q (R; w)-Q 00 =( ^ttR 3 ) Aw (p ) + (4^) - (p ° P °° )2 
+ (AttR) wl (2Aw (p ) - w - 2wi) + O (R°) 



(B5) 



+ 60w 2 ) 



(B6) 



Of course, we could solve the equation for the width exactly since it is simply a fourth order 
polynomial in w, 



= -K (p - Poo) 2 + 
Aw (p ) - 2wi 



K (p - p c 



+ 



R 



w 



3R 2 

6W2 — 2Aw (po 
R 2 



+ 2w — 2Aw (po 



w 



(B7) 



We can also express the free energy in terms of the equimolar radius 

Re = (± {R 4 -(R- -) 4 )) = R- 1 -w + ^(l- 6 Ql ) ^ 
+ ^(4a 1 -12a 2 + l)(^) 2 + 0(^ 



(B8) 



"Wo 

R 



we find that 



D r> > w <> > w o /a i\ w ° , ^o(6a 2 -5ai) ^ w 
i? = ^ + T + 12 (6ai - 1) ^ + 12 



R f 



O 



w 
Re 



(B9) 
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giving 

\2 



(ife; «,) - flee = Auj (p ) + U* f * ( *~ 2 Poo) + J (4^1) (BIO) 



+ u,2 ( 6 ^(Po 2 Poo) 2 + (13 + 3fli) Aw _ g _ _ 12 _ \ 2tt^ + Q ^ 
u> / 3 
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